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REACTIVITY  COEFFICIENT  MEASUREMENT  OF  BUCKLING 

Kenneth  W.  Downes 
Herbert  J.  Kouts 

Introductions  The  current  series  of  exponential  measurements  on  light  water 
moderated,  slightly  enriched  uranium  rod  lattices  is  intended  to  provide  typical 
reactor  parameters  for  such  assemblies  in  the  clean,  cold  state.  So  far,  ex¬ 
trapolating  to  operational  conditions  in  a  going  pover  reactor  (with  fission 
production  buildup,  high  temperature,  structural  poisons,  partial  fuel  burnout, 
Pu  buildup,  etc.)  must  be  done  by  calculations  which  are  first  made  to  predict 
the  clean,  cold  results  correctly.  Of  course,  it  is  desirable  to  make  such 
estimates  of  the  predicted  behavior  as  sound  as  possible* 

In  order  to  provide  some  experimental  basis  for  predicting  the  effect  of  25 
burnup  and  49  buildup,  we  are  planning  to  supplement  our  present  work  with 
measurements  in  assemblies  fueled  partly  by  Pu*  We  intend  to  measure  as  many 
quantities  as  possible,  hoping  to  get  critical  masses,  and  to  find  out  pertinent 
details  of  the  neutron  economy. 

We  have  also  scheduled  some  measurements  of  a  similar  nature  in  assemblies 
fueled  by  mixtures  of  25  and  thorium.  These  should  shed  light  on  the  utility 
of  such  &  system  as  a  25-23  converter. 

Since  fabrication  of  such  special  rods  is  expensive,  we  have  been  trying  to 
find  ways  to  measure  the  quantities  we  are  interested  in,  using  only  small 
amounts  of  fuel.  This  report  gives  the  results  of  a  measurement  of  buckling 
done  by  finding  the  reactivity  coefficient  of  an  unknown  fuel  in  terms  of  a  known 


Method:  Basically  the  measurement  consists  in  replacing  the  center  of  a  repro¬ 
ducing  assembly  by  a  small  reproducing  region  haying  an  unknown  buckling.  The 
critical  mass  is  measured  before  and  after  the  substitution,  and  the  change  in 
critical  mass  is  used  as  a  measure  of  the  unknown  buckling.  The  precise  value 
of  the  buckling  could.  In  principle,  be  found  in  two  ways.  The  measurement 
could  be  calibrated  by  using  in  place  of  the  unknown  a  series  of  fuels  of  diff¬ 
erent  but  known  B2;  the  unknown  value  could  then  be  found  from  a  curve  of  B2 
versus  change  in  critical  mass.  If  the  unknown  center  region  is  large  enough, 
it  may,  on  the  other  hand,  be  possible  to  determine  the  buckling  from  a  two- 
region  calculation. 

In  the  test  measurement  ve  have  made,  the  known  region  contained  .6 *  diameter 
uranium  rods  enriched  to  1*155?  25,  and  the  "unknown"  fuel  consisted  of  19  *6" 
diameter  rods  of  uranium  enriched  to  1.02(75?.  The  moderator  ws  light  water, 
in  effectively  infinite  water  reflector  was  used.  The  bucklings  of  both  enrich¬ 
ments  at  the  water-to-metal  volume  ratio  used  (3:l)  were  known  beforehand,  and 
the  purpose  of  the  measurement  was  to  find  out  if  siarole  two-region  analysis 
would  be  adequate  for  predicting  the  of  the  1.027?  enrichment  rod  region. 

The  critical  masses  for  thp  two  cases  were  found  by  a  method  described  in  BNL 
Log  No.  C-7592.  The  uranium  was  loaded,  with  suitable  precautions,  in  the  pre¬ 
sence  of  a  neutron  source  until  the  multiplication  of  neutrons  from  spontaneous 
fissions  was  great  enough  to  permit  flux  level  measurements. 

At  this  point  the  neutron  source  was  removed  and  the  neutron  flux  in  the  assembly 
was  measured  as  rods  were  removed  from  the  periphery.  The  result  of  the  measure¬ 
ment  was  then  an  experimental  plot  of  the  local  neutron  flux  as  the  fmiction  ef 
the  loading. 
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Because  of  spontaneous  fissions  in  a  simple  one-region  assembly  the  flux  density 
in  the  core  has  the  form  (we  include  the  effect  of  the  reflector  through  an 
assumed  constant  reflector  savings): 


0  =  2  c  sin  JLlZ.  J  L 
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Q  =  source  of  spontaneous  fission  neutrons 
h  =  height  of  the  assesnbly 
5  k  =  successive  zeroes  of  J  (u) 

T^t  -  transport  rafp  ° 

A  =  reflector  savings 

For  the  nearly  critical  lattices  ve  have  used,  only  the  fundamental  mode  j,k  =  1 
is  important,*  so  that 
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(R  the  critical  loaded  radius)  for  a  critical  finite  cylinder.  Thus 
c 

constant  •  sin  jlj.  jo  (SjlL) 

0  -  h  \K  +  A/ 


-  ll 

iR  +  a/  \A+ A 


The  thermal  flux  is  measured  at  a  fixed  point  on  the  central  axis.  Thus  a  plot 
of  1/0  T8.  /  1  r  is  a  straight  line  .intersecting  the  axis  at 

ivN 


(R  ft  +X)2 

c 

1.  The  fluxes  measured  had  .01$  harmonic  content. 
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Of  course,  the  procedure  is  essentially  that  of  the  critical  assembly.  But 
since  the  lattices  are  not  loaded  all  the  way,  the  analytical  justification 
for  extrapolating  the  flux  plots  to  critical  must  be  made. 

In  the  two  region  assembly  the  experimental  procedure  was  precisely  the  same. 
Analysis  of  the  experimental  results,  however,  should  strictly  be  done  by  a 
two-region  calculation.  We  did  not  do  this.  What  we  did  was  to  establish  that 


the  experimental  riot  of 


is  in  this  ease  a  straight  line 


also,  and  assume  that  the  intersection  of  this  straight  line  with  the  axis 
correctly  indicates  the  critical  mass. 

Experimental  procedure;  The  lattices  were  loaded  up  to  a  suitable  with 

appropriate  flux  monitors,  safety  rod,  and  trip  circuits,  and  with  flux  levels 
maintained  by  multiplication  of  neutrons  from  a  Po-Be  source.  At  this  point 
the  source  was  removed  from  the  lattice,  and  taken  away  from  the  vicinity.  A 
small  counter  was  inserted  in  the  core,  and  Count  rates  were  recorded  for 
several  smaller  loaded  radii. 


The  maximum  count  rates  from  spontaneous  fission  multiplication  were  about  250 
c/m  and  a  range  of  about  a  factor  of  3  in  neutron  multiplication  was  used. 

Analysis:  As  has  been  mentioned  earlier,  all  loadings  were  done  with  an  effect¬ 
ively  infinite  water  reflector.  The  simplest  calculation  of  the  critical  mass 
of  the  assembly  with  the  center  replaced  with  an  "unknown*  is  done  with  ordin¬ 
ary  one-group,  two-region  diffusion  theory,  with  the  effect  of  the  reflector 
being  lumped  into  an  assumed  constant  reflector-savings. 

In  this  case,  if  ve  let  the  subscripts  1  and  2  represent  the  inner  and  outer 
region,  respectively,  the  pile  equations  are 


V  ^  =  0 
y2  4 +  E2 1*2  - 0 
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The  solutions  are 


Bf  =  a  J  (ur)  +  b  T  (ur)  I  cos  TT  a 
1  L  o  o  J  h 

fS2  =  j^c  JQ  (vr)  +  d  Tq  (vr)J  cos  -X_B 


=  u2  +  ^  TT  j 
B2  =  v2  +  (  IT  \ 

2  'sTl 


The  boundary  conditions  are 


^  =  finite  at  r  =  0 
^  =  02  at  r  =  %  (interface  flux) 

-  ^2 _ at  r  =  (interface  current) 

dr  dr 


$2  =  0  at  r  =  Rg  (extrapolated  boundary) 

In  the  third  boundary  condition  ve  have  assumed  the  equality  of  the  diffusion 
coefficients  for  the  two  regions.  In  general,  of  course,  we  may  not  assume 
this,  but  the  two-region  assembly  we  are  analyzing  here  very  closely  satisfies 
this  condition. 

Insertion  of  the  solutions  (8)  into  the  boundary  conditions  (10)  yields  readily 
the  secular  equation 

u  J  (uR.)  =  vfj,  (tR1)Io(tR?).I1(tR1)J0(TR2)| 


J0  (u  %) 


Tq  (  ▼  R2  ^  Jo  ^  v  ®1  ^  •  To  (  v  R1  ^  Jo  ^  v  R2  ^ 


Supposing  that  R1  and  R2  are  town®  and  that  ▼  is  known  from  (9) ,  we  may  solve 
for  u.  Then  use  of  (9)  gives  B^,  the  unknown  buckling. 
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Ve  hare  also  carried  through  the  two -group,  two-region  analysis.  For  the  ex¬ 
perimental  case  ve  are  analyzing,  the  tvo-region  diffusion  coefficients  and 
diffusion  lengths  are  very  nearly  the  same  in  the  outer  and  inner  region.  If 
we  assume  they  are  exactly  equal,  then  the  two-group  secular  equation  reduces 
again  to  (10) . 


The  experiment  and  its  analysis  are,  stepwise,  as  follows: 

(1)  A  valtie  of  Ms  assumed  from  previous  exponential  experiments. 

(2)  The  neutron  flux  resulting  from  multiplication  of  the  spontaneous 
fission  neutrons  is  measured  for  the  one-region  assembly  as  a  function 
of  loaded  radius  R. 


(3)  A  straight  line  plot  of  1  f$  vs 


is  made,  and  the  value  of  * 


the  critical  radius  Rc  is  found  from  the  intercept  of  this  straight 
line. 

(4)  is  evaluated  as 


(5)  The  inner  region  of  this  assembly  with  radius  R. ,  is  then  replaced 
by  the  unknown,  and  steps  (2)  and  (3)  are  repeated.  The  critical 
radius  in  this  case  is  R2. 

(6)  These  values  of  B^,  1^,  R2  are  then  used  in  equations  (9)  and  (ll) 
to  give  P|. 


Results :  The  value  of  a  used  was  6.6&  cm.  This  value  was  found  from  previous 
exponential  experiments  with  the  3sl  volume  ratio  lattice  of  1.15%  enriched 
rods  in  ordinary  water.  The  experimental  inverse  flux  plots  of  the  spontaneous 
neutron  multiplication  are  shown  in  figure  1.  It  is  evident  that  the  points  in 
both  cases  are  well  represented  by  straight  lines.  This  is,  of  course,  to  be 
expected  for  the  one-region  assembly,  because  equation  (5)  predicts  this  behavior. 


-  6  - 


A  straight-line  fit  to  the  two-region  data  seems  equally  justified  by  the 
appearance  of  the  plotted  points* 

The  critical  radii  for  the  two  cases  are  found  from  the  intercepts  of  least- 

squares  fit  straight  lines*  They  are 

R  =  37.498  ±  o007  cm  (one-region) 
c 

R2  =  38.088  ±  .06 L  cm  (two-region) 

with  the  error  limits  established  by  the  standard  deviations  from  the  best 
straight  lines.  Actually,  the  accuracy  of  the  measurements  is  probably  closer 
to  the  latter  value}  the  points  for  the  one-region  measurement  fit  the  straight 
line  better  than  would  he  expected  just  from  the  ±  1%  counting  accuracy. 

We  find 

v  =  £t4P.4.8 —  =  6.4131  x  10"2  cm“l 

*c 

The  chosen  value  of  %  is  based  on  the  volume  ratio  (3  si)  and  the  number  of 
rods  used  in  the  central  region  (19).  We  get 

R^  =  6.8325  cm 

The  value  of  u  is  obtained  from  inserting  these  values  for  R^  Rj,  and  v  in  (9) 
and  (11).  It  is 

u  =  5.2755  x  10“2  cm*1 

The  height  h  of  the  assembly  is  taken  from  geometry  to  be  133-6  cm.  Thus  we 
get 

B2  =  (5.2755  sc  10-2)2 

=  33.36  ±  ,50  x  ltr*  cm*2 

with  the  error  estimate  based  eventually  on  the  least  squares  fits  to  the  in¬ 
verse  flux  plats.  This  value  is  to  be  compared  with  that  determined  from  pre¬ 
vious  exponential  experiments: 
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The  two  values  agree  to  and  are  within  the  nutual  error  limits. 

Certainly  a  large  part  of  this  excellent  agreement  can  he  associate-4  with  the 
similar  neutron  slowing-doun  and  diffusion  characteristics  of  the  two  regions. 
If  they  were  not  similar,  the  two-group  analysis  would  not  give  the  same  value 
of  as  the  one  group  analysis,  and  results  would  be  more  uncertain.  We  may 
expect  then  that  our  scheduled  measurements  with  alloyed  \%  Pu,  99%  U  rods 
will  produce  good  values  of  B^.  Those  made  with  the  thorium-uranium  rods  may 
not  be  as  trustworthy.  We  expect  to  find  out  something  about  their  validity 
by  measuring  with  two  enrichments  of  uranium  in  the  outer  region. 
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